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^\1 . We present a general linear response description of membrane adhesion at rough or chemically 

structured surfaces. Our method accounts for non-local Van der Waals effects and contains the more 

approximate (and local) Deryagin approach in a simple limit. Specializing to supported membranes 

ry\ ■ we consider the effects of substrate structure on the membrane adhesion energy and conhguration. 

Adhesion is usually less favorable for rough substrates and the membrane shape tends to follow 
that of the surface contours. Chemical patterning, however, favors adhesion with the membrane 
configuration being out of phase with the surface structure. Finally, considering a surface indented 
with 'V'-shaped trenches, we show that our approach is in fair agreement with an exact numerical 

fi^ , solution. 
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I. INTRODUCTION 



g : Supported membranes strongly adhere to substrates and lie typically at separation distances of between 10 A and 

fj ' 40A |l|. Such small values have lead to their adoption by the biotechnology industry \^ and, in particular, given 
them an important role in the development of biosensors. Supported membranes enable one to biofunctionalize an 
inorganic surface and can provide an ultra-thin, highly electrically resistant layer on top of a conducting substrate. 
They provide a means of immobilizing proteins with a well-defined orientation and in a non-denaturing environment 
.^ Q . If these proteins are receptors then one can use electrical or optical means to detect or "sense" the binding of 

pg ligands to the receptors [Q . 

OQ Supported membranes can be formed by the spreading of a bilayer over a substrate, vesicle fusion taking place at a 

On ' substrate or by lipid monolayer transfer using a Langmuir-Blodgett technique H]. However, in nearly all applications 

>~j ' the substrates used are not simply planar and homogeneous but are patterned, either chemically HJq| or geometrically 

|g|,0. The theory of membrane adhesion has typically concentrated on adhesion at ideal planar surfaces ^. In a 

recent paper [p|, we have provided a simple description of the adhesive properties of membranes at rough surfaces. In 

ti . this article, we would like to present a more general approach to membrane adhesion which includes the possibility 

j^ ' of chemical patterning. 

We begin with a summary of the basic assunrptions of our model and of the intermolecular interactions involved. In 



'"^ ■ Sec. [I], a planar, homogeneous substrate is considered as the starting point for the linear response theory that follows. 
^ We show how the simpler approach of Ref. [|| is incl uded in the present work and then proceed to consider several 
O illustrative examples in Sec. |Vl[ Finally in Sec. VII, our analytical description is compared and contrasted with a 



complete numerical solution. Throughout, we emphasize the effect of substrate roughness and chemical heterogeneities 
^ on the adhesive properties of the supported membrane. 

X: 

C^ ■ II. THE FREE ENERGY 

To begin we consider a membrane supported at a substrate that can be either geometrically structured (non- planar). 
Fig. 0(a), or chemically structured (patterned with different chemical compounds). Fig. n|(b) and (c). For example, 
a surface can be chemically structured by depositing different chemical layers, see Fig. 1Kb), or adjoining different 
chemical surfaces together to make a columnar structure. Fig. |l|(c). For a membrane adhering to such a surface this 
patterning will greatly influence the membrane configuration and adhesion energy. Inspired by advances in the theory 
of wetting [|lO| , we adopt a general mean-field approach in which the configuration taken up by the membrane is one 
that minimizes the free energy. In order to find this optimunr configuration we first discuss the form of the free energy 
functional. 



If the membrane has an elastic modulus k and tension cr, then its bending energy can be described by the functional 
ill 



dSy/g 



a+-n{2Hy 



(2.1) 



where the integral is over the membrane surface, g is the determinant of the metric, H — [ci -\- C2)/2 the mean 
curvature, and ci and C2 the two principal curvatures. We have ignored here the Gaussian curvature contribution as 
only a membrane with a fixed topology (flat on large lengthscales and of infinite size) is considered |12[| . Throughout 
the paper we choose to work in the Monge representation. Letting p = (x, y) be a two dimensional planar vector, the 
heights of the surface and membrane above some reference p-plane are Zs{p) and h{p), respectively (see Fig. |l|). 

To account for the interaction of the membrane with the substrate, we include a potential term, V{h] Zs^e), in 
the free energy. As already mentioned, Zs{p) accounts for the substrate's geometrical structure and describes its 
surface configuration, while e(r) denotes any chemical inhomogeneities. The potential can have a number of different 
components [HSUmI. For our case, the most important of these is the Van der Waals contribution, which is given by 



Kdw(/i;2s,e) 



W(h-z,,e)-W(h + 5-z,,e) 



(2.2) 



where 5 is the membrane thickness and is typically around 40 A. Due to the bilayer nature of the membrane, two 
terms involving the Van der Waals potential, W{h] Zg, e), are necessary; in particular for supported membranes where 

For a thin fluid film of thickness h{p) resting on an inhomogeneous solid, one can sum over all possible pair 
interactions between the molecules in the upper half space, capped from below by the surface z = h{p), and those in 
the lower half space, capped from above by z = Zs{p), to show that W{h; Zg) satisfies [Q 



Wih;zg,e) = I dz j d^ p' 
Jh(p) 



^s(p') 



dz' wo{r-r') 1 + e(r') 



with 



Mr) = -J -6 



(2.3) 



(2.4) 



The latter models non-retarded Van der Waals interactions. Equation (2.3) contains a position dependent Hamaker 
constant 



^(r) =Ao 1 + e(r) 



with Aq the average value 



Ao 



and e the (small) deviation around this average, 

.(r) 
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(2.5) 
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We emphasize that W{h; Zs,e) is a functional of both Zgip) and e(r). 

If h{p) is set to a constant value, say ho, and both Zg and e vanish, then ( |2.3| ) becomes the familiar 

Wiho;0,0) = Wo{ho)^4r-T2 



(2.8) 



(see Ref. jlj]), which is just the Van der Waals potential between two planar semi- infinite bodies held a distance ho 

apart. 

Equation ( |2.3| ) provides an attractive interaction and, for the case of a supported membrane, this is chiefly balanced 
by hydration forces. The hydration potential obeys 



Vhyd{h;zs 



be- 



-a{h-Zs) 



(2.9) 



where b has units of surface tension and a is an inverse length of typical size a~^ ~ 2-3 A. Due to the very short 
range nature of the hydration interaction, we include the dependence on the substrate structure with a simple local 
approximation and so T4iyd is just a function of the local height h(p) — Zs{p). The origin of hydration forces is still 
under debate (l4| but they are generally believed to have some steric contribution. Consequently, while b in general 
is position dependent we believe that this is a relatively minor effect and so choose to keep the simple form of (2.9). 
The total potential is then 



V{h] Zs, e) = Kdw(/i; 2^8, e) + Vhyd(/i; 2s) 



(2.10) 



though one could consider more complicated scenarios involving, for example, electrostatic forces. Summing all these 
contributions we can write the total free energy as 



T[h]= f d^pWi + {vhy 




Vh 



VT+WW 



+ V{h;zs,e) 



(2.11) 



where we have explicitly written out the curvature and tension terms in the Monge representation. 

One of the most relevant quantities in experiments is the membrane adhesion energy. Within our general mean-field 
approach, the optimal height of the membrane is that which minimizes (2.11) and the value of the free energy when 
the membrane takes up this optimum configuration, J^min, leads to a natural definition of the adhesion energy per 
unit area 



U=- 



T„ 



(2.12) 



Here, S'o is the total area of the projected reference p-plane, Sq — J dP p and we have subtracted off a membrane 
tension term. Doing so conveniently defines the adhesion energy so that a completely flat membrane does not have a 
tension dependent contribution; for a membrane infinitely far from the surface U will then vanish. Notice that (2.12) 
implies that an attractive surface will have a positive adhesion energy. 



III. PLANAR AND HOMOGENEOUS SUBSTRATES 

Our results are obtained by analytically expanding the free energy around its value taken for a planar, chemically 
homogeneous substrate. Therefore, we briefly review the results for such an ideal surface. 
For this case, e = Zg = 0, and the Van der Waals interaction (^.3|) simplifies to (^.8[), and so 



Kdw(/i; 0, 0) = - Wo{h) -Wo{h + 5) 



(3.1) 



where throughout we use the subscript zero to denote adhesion at both chemically homogeneous and flat substrates. 
Here, a is a fundamental lengthscale in our problem 



\2Tia 



and is provided by the ratio of the Hamaker constant, see (2.8), and the membrane tension p5[ 
We find that the membrane adopts a flat configuration, h[p) = /iq, which obeys 



dV_ 
'dh 



= K'dw(^o) -abe 



-aha 



(3.2) 



(3.3) 



from (|2J). 

The adhesion energy in this case is simply given as the negative of the interaction potential. From (2.11), J-^in = 
So [<J + V{ho; 0, 0)] and so ( ^.12| ) implies that 



Uo = -VihQ;0,0) 

= -[Kdw(/io) + 6e-"''°] 



(3.4) 



By definition, Uq is positive for all sufficiently attractive potentials, V. Equations (3.3) and (3.4) provide the funda- 
mental quantities upon which our perturbation theory will be built. 

In order to allow (semi-quantitative) comparison with experiment and to give some idea of the magnitude of the 
quantities involved, we would now like to specialize to a particular choice of our model parameters (we opt again 
for those chosen in Ref. B), see Table |[ Typical experimental values of a and n are 1.7 x 10~^ Jm~^ and 35T, 
respectively [|l^. We set the Boltzmann constant to unity and so at room temperature T = 4.1 x 10^^^ J. Choosing, 
Ao = 2.6 X 10^21 jjj^-2 ^ 0.63 r ||l|l, implies that the Icngthscale a ~ 49.3 A and, from (|!|), ho ~ 0.61a ~ 30 A in 
agreement with measured values using specular reflection of neutrons [ |l7| . The two parameters used here to specify 
the hydration force, see ( |2.9| ), are 

6 = 0.93Jm-2 ; a~i = 2.2A (3.5) 

which are in accordance with those measured in Ref. pq | . 

The potential experienced by the membrane, ( ^.10 ), is sketched in Fig. |[ From (3.4), one can see that 



Uo 



0.298cr = 5.07 X 10"^ Jm"^ 



(3.6) 



At this point, it is also worth discussing the other lengthscales which will appear in our treatment. Defining v as 
the second derivative of the potential calculated at the minimum, h = Hq, 



^Vih;0,0) 



several correlation lengths can be extracted 



together with 






il = <^h ; i\ = i^h 



e' = «/^-dc 



(3.7) 



(3.8) 



(3.9) 



which describes the crossover between the tension and the rigidity dominated regimes. Their values for the experi- 
mental scenario described above are given in Table E (which also lists all the other model parameters). 



IV. LINEAR RESPONSE THEORY 



To carry out a perturbation theory for rough and heterogeneous substrates, we assume that Zg and h — ho are small, 
i.e. Zs '^ h — ho <S^ ho^ and that any of their products and derivatives are also small. A (functional) Taylor expansion 
is then performed which is a fairly standard, if long, calculation. 

To simplify our presentation and ease the algebra, we will assume that the chemical structure is such that e can be 
factorized, i.e. 



e(r) ^(f>{p)g{z) 



(4.1) 



for some functions and g and the substrate surface is given by Zs = 0. Such a factorization while including the 
layered (constant (p) and columnar (constant g) structures shown in Fig. n^ does prevent us from considering surfaces 
which are both rough and chemically inhomogeneous. Consequently from this point on, we will specialize to either 
rough or chemically patterned substrates. 



A few more definitions are in order; first of all, we notice that w, given by (3.7), can also be shown to obey 

V = V{:^^{ho) - j d?p'[wo{p\ ho) ~ woif', ho + (5)} (4.2) 



using (2^) and (2^). The kernel functions (this choice of nomenclature will become clear later) 



K{p) = -- [wo{p,ho) ~wo{p,ho + S)] 

V 



(4.3) 



and 



G{p) 



ha+S 



dzg{ho - z)wo{p,z) 



will also prove useful. 

Then, expanding J-[h] in ( ^.11 ) to second order and taking advantage of (0.3|), we find 



l^^l^o+^J d'pUi'^hf + e^^^hf + ih- ho? 



2{h~ho) Vzs+ / d^p'K{p')z,{p + p') +z^{p) 



for rough surfaces and 



l^^l^^+^j d'pUii^hf + eA^^hf +{h- hof 



2{h-ho) / d^p'G{p')cf>{p + p') 



for chemical structure, where J-'q — J-^[h()\ is the /i-independent term in the expansion. 



To^ <J 
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In the case of chemical patterning, Tq contains an additional term 

which can be made to vanish by choosing the /?-plane such that 

d^p^{p)=0 

For rough substrates the p-plane is chosen so that (zs) = / d'^pzs = 0. Here, V is given by 



V 



V^Aho) 



(4.4) 



(4.5) 



(4.6) 



(4.7) 



(4.8) 



(4.9) 



(4.10) 



To find the optimum profile, we need to minimize (4.5) or (4.6) with respect to h{p), i.e. SJ-/6h — 0. The resulting 
Euler-Lagrange equation is non-local in p but linear in h, 



^y'~ev' + i 



i^h{p) - ho 



Jd'p' 
Jd'p' 



K{p' - p)zs{p') + Vzs{p) rough 

G{p' — p)(f>{p') chemical 



(4.11) 



where the role that K and G play as kernel functions becomes clear. Equation (4.11) is the starting point for our 
linear response profiles. It is a fourth order non-homogeneous linear differential equation where the heterogeneity of 
the substrate enters in the non-homogeneous term. Due to its linear nature, the solution can be written down in 
Fourier space. Defining for any function f{p) its Fourier transform 

(4.12) 



we find, via the convolution theorem, that Sh{p) =^ h{p) — ho obeys 



K(q)+V 



s(q) 



dh{q) = < 



H-e<?"+?J'7'' 



G(q)0(q) 
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(4.13) 



The Fourier transform -fir(q) of the kernel function ii'(p) can be calculated using the result 

e^'^P . / g \" i^™(g/i) 



d^p 



(p2 + /i2)m+ 



r-'© 



r(TO + i) 



(4.14) 



where Km{x) is the modified Bessel function of the second kind of order m. Non-retarded Van der Waals interactions 
are obtained by setting to = 2 in the above equation (from (2.4)). Then the kernel becomes 



K{q) = K{q) 



2„2„2 



e.q'a 



K2{qho) K2{qhQ + qS) 



''o 



{ho + sy 



(4.15) 



Notice that in the limit of q tending to zero, limc^^o Q^K2{q) ~ 2, which implies that 



i^(o) = -e>^ TT 



ht {ho + S)^ 



1-V 



(4.16) 



from i^) and ( pO| ). 
Similarly, G'(q) satisfies 



which then implies that 



G(q) 



t2„2„2 /-/io+t) ^/u ^\ 



G(0) = -Ca' 



'"^' .. 9{h. - z) 



dz 



For the case when g{z) = 1 (a columnar solid) 



G(0) 



ea' f 1 



hi {ho + 5f 



<0 



(4.17) 



(4.18) 



(4.19) 



which will be needed in Sec. VI 



We can also calculate the adhesion energy in Fourier space by Fourier transforming (4.5). Using the solution of the 
Euler-Lagrange equation, ( 4.13| ), and the definition (2.12) then gives 



U 
a 



Uo 



k/721^ l^~s(q)p- 



(^K(q)+vJ£,(q) 



i+qHl+qH 



|G(q)0(q)r 
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(4.20) 



V. THE DERYAGIN APPROXIMATION 



In a previous paper [^, we looked quite extensively at the adsorption of membranes on rough substrates and 
throughout used a Deryagin-like local approximation ||l3]. In this section we would like to show how this is included 
in our present, more general, linear response approach. 



The Van der Waals potential, (^.2|), due to its functional (non-local) dependence on the inhomogeneities, provides 
most of the difficulties in any analytical method. These complications are neatly circumvented by the Deryagin 
approximation. Equation (|2.3|) can, by a simple change of variable, be re-written as 



Wih;zs,e) 



j;:^^^dzJd^p'f^^dz'wo{p',z^z')[l + ^{p + p')giz') 
where we have adopted ( |4.l| ). The Deryagin method amounts to replacing 



rough 
chemical 



Zs{p + p') ~* Zs{p) 
Cj){p + p') -. CJ3{p) 



rough 
chemical 



(5.1) 



(5.2) 



which removes the func tional character of (5T) and so neglects almost all non-local effects (we still retain the integral 
over g{z')). Once ( |5.2| ) has been performed, one can (Taylor) expand the free energy in powers oi h — ho and Zg as 
before. 



However, such an approach turns out to be equivalent to replacing the kernel functions, (O) and (4^), in the linear 
response theory by Dirac delta functions. 



K{p^p')^K{0)S{p-p') 
G{p-p')^G{0)S{p-p') 



(5.3) 



which transparently shows the local character of the Deryagin technique. Equation ( p.3[ ) implies that the Fourier 
transforms of the kernel functions are now sim ply co nstants. 
Consequently, the Euler-Lagrange equation, ( 4.11 ), becomes 



(dv^ - ev' + i)sh{p) = 



Zs{p) 

G(O)0(p) 



rough 
chemical 



where we have used (4.16). In Fourier space the solution is 

I 1+^ 
Shiq) 
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i+iiq^+iU' 
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(5.4) 



(5.5) 



with an adhesion energy given by ( 4.2Cl|) with the substitution (5.3) 

Uo. i_ r d'^q. 

a 2So J (27rf„ 



.2t2 I „4t4 



u _ 

(7 



Ug , 1 f d 
a "'" 25o J (27r 



d^q 



.r 



i+q'^e+qH 



|G(o)0(q)r 
i+q^-e+q'^ii 



rough 



chemical 



(5.6) 



For the case of roughness, these are exactly the results obtained in our previous paper Q. We note that there is a 
difference in sign between the two cases; while 6U = U ~ Uq < for the ro ugh ca se, 5U > for chemically patterned 
surfaces. This is one of our main observations and will be discussed in Sec. VI C . 



VI. PREDICTIONS OF THE LINEAR RESPONSE METHOD 



The non-local perturbation method is embodied by (4.11) and (4.20). In this section we compare and contrast this 
approach with the simpler local Deryagin approximation. As both methods are based on perturbation theory, we 
have the caveat that 



U/Uo 



1 



(6.1) 



for the method to be valid, i.e. that the perturbation correction is much smaller than the term it improves on. The 
fundamental difference between the approaches is in the treatment of the Van der Waals potential; in the Deryagin 
approximation there is only a p urely local attraction while the linear response includes non-local collective effects. 
The Euler-Lagrange equations, (4.11) and (5.4), most clearly illustrate this. 



A. Chemically structured surfaces 



For the types of structure considered, ( |4.5| ) gives the general linear response free energy. One can define a bending 
(BE) and potential energy (PE) contribution to this by simply letting the bending energy be that part which vanishes 



when a and k are set to zero and the potential energy that part that remains. Then, using the solution, (4.f3), we 
find 

PEeW = -/|^(i±^|G(q)^(q)P (6.2) 



r|G(q)0(q)P (6.3) 



and 

dPcj Q 

(27r)2 (1 + Q): 

with Q — ^^q^ + ^t.^"^ (^ constant contribution, y(/io;0,0), has been ignored in the potential energy term). Here 
and throughout the rest of the paper, we fix g{z) = f for clarity and so consider only columnar solids. The reader 
interested in layered solids should consult Ref. ||T^] whose results can be simply extended to membranes. 
Using the fact that 

> G[q) > G(0) (6.4) 

for all g, and remembering that the Deryagin approximation is recovered when the kernel function, G{q), is set to 
its value at q = 0, one can see almost by inspection that the membrane potential energy will be higher than in the 
Deryagin case while the bending energy will be lower. In fact, the increment in the potential energy will be greater 
than the bending energy decrement and we can therefore expect the surface to appear less attractive, due to non-local 
effects, with a correspondingly lower U. 

A sinusoidally patterned surface, translationally invariant in the y-direction, is the simplest choice with which to 
illustrate this behavior. Setting 

4>{p) = Ac sin(pca:) (6.5) 

with Ac being the amplitude of the Hamaker coefficient oscillations and 2t:/pc their period, ( 4.13| ) then implies 



AcG'(pc) sm{pcx) 



h{p) = /lo + ^\ ,2'2_ufLi (6-6) 



where G{pc) < from (6.4). The adhesion energy is also easily calculated 

U^U^l_ Ag|G(Pc)P 



(6.7) 



In Fig. 0(a), this adhesion energy is plotted against wavenumber pc- One can see that there is good agreement 
between the Deryagin and linear response approaches. As the period of the sinusoidal variation increases the membrane 
is less and less able to respond to the structural variations and for high pc the adhesion energy takes its planar value 
once more. Any chemical structure has been effectively washed away. 

B. Rough surfaces 

In this case, the two contributions to the adhesion energy are 

i i^rough - y (27r)2 1^'''-^-' 
and 



PErough = J ^ \zs{q)\' |1 - |i^(q) + ^\ \^^Q^2 > (6-8) 



BErough = / 7^ — %^|(i^(q) + V)f.(q)p (6.9) 



(27r)2 (1 + Q)2' 



As 



Kiq)+V>K{0) + V 



(6.10) 



for all q, the complete opposite behavior results with the bending energy increased by the non-local contributions and 
the potential energy decreased below the decrement to the bending energy. The surface becomes more attractive and 
the adhesion energy inceases above the Deryagin result. 

Looking at a chemically homogeneous but geometrically corrugated surface we choose 



Zsip) = As sin{psx) 
so that the surface corrugations have an amplitude of As and a period of 27r/ps. Equation (4.15|) gives 



with K{ps) > from ( 6.10| ), while U obeys 



(6.11) 



(6.12) 



(6.13) 



This latter result is plotted in Fig. g(b) and deviates substantially from the Deryagin prediction as Ps increases. 
Non-local effects are important and can strongly decrease the membrane potential energy. In particular, notice that 
for small Ps the adhesion energy is i ncre ased above the value taken for the planar situation. This can never occur 
in the Deryagin approximation (see (^) and Ref. ^). The additional Van der Waals contribution accounted for 
substantially raises the membrane potential energy. For large ps, positive non-local effects "saturate" i.e. K{ps) 
plateaus, and U starts to decrease for greater p values, see ( |6.13 ). In general, the membrane finds it more difficult to 
adhere to the rough surface and the adhesion energy will asymptotically (large Ps) tend to a constant value less than 
unity (though for the high value of As chosen our perturbation theory is not sufficiently accurate to capture this). 
The membrane potential energy eventually becomes positive while its bending energy vani shes as Ps ^ cx). 

Again, it is important to reiterate that the our approach is strictly only valid when (6T) holds. 
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C. Comparison between chemically structured and rough surfaces 



Consider once more (|4.13D , with Q = ^"^q^ + Q.q 



5/i(q) 



K(c^)+V i-.(q) 



1+Q 
G(q)0(q) 
i+Q 



rough 
chemical 



(6.14) 



Two points are worth making; (i) from (6^) and ( 6.10| ). It can be seen that the membrane amplitude is increased and 
decreased by non-local effects for roughness and chemical structure, respectively. This is reflected in a corresponding 
change in the bending energy (compared to the Deryagin results), (ii) That due to the different signs oi K + V and 
G, see (6.4) and ( 6.10| ) again, the membrane profile is always in phase with the surface contour of a rough substrate 
but is exactly out of phase with surface structure arising due to chemical variation. A membrane adhering to a rough 
surface will, in order to maximize its potential energy, try and follow that surface as best it can (limited only by 
the resulting bending energy cost). Similarly, the membrane will follow the hills and valleys of a substrate potential 
generated by chemical patterning. In this case, however, a local increase/decrease in the Hamaker constant, A, makes 
that region of the surface more/less attractive and so shift the membrane in/out. Note that the membrane would be 
in phase with the substrate for a repulsive Van der Waals potential [ pO[ . 

By using the Deryagin approximation, simple arguments describing an effective potential energy highlight the 
different adhesive properties of the two surfaces. In Fig. 0, a surface (translationally invariant in the y-direction) 
patterned with a square wave profile and a similar on/off chemical patchwork is schematically shown. For the former 
case, a definition for an effective potential energy is given by 



V,s 



JdxV{h-Zs) 
Jdx 



(6.15) 



which just equals V^ho) for the planar scenario, Zg — 0. When the substrate is square-corrugated this approximately 
becomes 

Vcs « I {Viha + Sh-c) + V{ho -Sh + c)} 

- V{ho) + liSh - cfV"{ho) + 0({5h - cf 

> V{ho) (6.16) 

as V"{ho) is positive {ho being a minimum of V). Thus even in this crude argument, one can see that surface roughness 
acts (at least if the roughness does not get too large when our perturbatio n ex pansion brea ks do wn) to increase the 
membrane potential energy — a result verified by the Deryagin version of (6^). Equation ( 6.16| ) shows that surface 



roughness is reminiscent of Gaussian thermal fluctuations and similarly to these acts to drive the membrane out of 
its potential minimum. 

When the wall is patterned with a periodic array of alternating chemical patches as in Fig. 0, one can use similar 
arguments to those given above to estimate the effective potential energy. Let the Hamaker constant A obey 

. _ J Ao(l + 0) dark patch , 

1 Aoil - (f>) light patch ^^'^'^ 



for constant and positive Aq and </), then the potential energy is, see ( p.5[ ) and ( 2.10 ), 

Vih)^V^dUh)-e{x) + V^yiih) (6.18) 

with e{x) — Izt (f> on the dark/light patches. Therefore 

JdxVjh) 

= ^[v^dwiho - 5h) ■ {I + (/)) + yhyd(/io - 5h) 
+Kdw(/io + Sh) ■ {I -(/)) + Vi,yd{ho + Sh)^ 
= V{ho) + ^6h^V"{ho) - ,p5hVU^{ho) + OiSh") (6.19) 

Notice here that the membrane's position is exactly out of phase with the surface structure. This leads to a negative 



contribution in (6.19) and one can see that the new membrane configuration, as the negative term is of order 5h and 
positive is of order 5h? , is likely to result in a net attractive contribution to the potential energy. Indeed, this can be 



verified by summing ( |6.2| ) and (6.3) which, as already mentioned, is always negative. 

In summary, a membrane generally adhers less favorably (relative to the planar and homogeneous surface) to a 
rough substrate and adopts a configuration that is in phase with the surface contours. A chemically structured 
substrate has a higher adhesion energy (more favorable) and leads to a membrane configuration exactly out of phase 
with the surface structure. 

VII. COMPARISON BETWEEN EXACT AND APPROXIMATE SOLUTIONS 

From these examples and those given in Ref. H, it is clear that the Deryagin approximation is certainly the most 
versatile if one wishes to obtain analytical results. However, it is also apparent that in some situations non-local 
effects can become important and in this section we compare the Deryagin result with an exact numerical solution. 

The scenario we choose to specialize to involves solely geometric inhomogeneities. We consider a chemically pure 
substrate made up of a regular array (in the x-direction) of 'V shaped trenches. These could be formed, for example, 
by etching silicon wafers |gl|]. In the other spatial dimension, i.e. the y-direction, the system is translationally 
invariant. For simplicity, the trenches are assumed to be symmetric about their lowest point and have a maximum 
width of d and a depth of Ad. See Fig. |^ for an example. 
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A. The Deryagin solution 



To find the Deryagin solution it is most easy to consider (5^) directly, which reduces to a one dimensional differential 
equation. This is 



with 



^'^-^'^ + ^'^'^(^) = "^(^) 



,{x)^2X Y^ <r^ 0{x + nd) - e{x + nd - d/2) 
+ {d - r^) 9{x + nd- d/2) - e{x + nd - d) 



(7.1) 



(7.2) 



where 9{x) is the Heaviside function and r^ = xmodd, i.e. the remainder of x on division by d. Assuming that the 
membrane has the same symmetry as the substrate, the boundary conditions are 



h'{nd) = h"'{nd) = 
h'ind + x) =h"'{nd + x) = 



(7.3) 



for all integer n. 

Due to the periodicity of Zs{x), we need only solve (7.1) for < x < d/2, where Zs{x) — 2Xx, and then reflect 
and/or translate this solution to obtain the full membrane configuration. We find in this region that for p = {x,y), 



2\ti^ Tl^ 
hip) — hQ + 2\x „ ^ ~ \il^+ix) — ib~{x)] 



with 



sinh 



■)p±{x) = 



{d/4 - x)r/± 



77^ cosh 



(77±d/4 



(7.4) 



(7.5) 



The ?7± come from the factorization of the operator in ([7T|), see Ref. 0, and are 



1 ± VI - m^n)' 



(7.6) 



where ^ is given by (p^). The height profile of the membrane, h{x) is plotted in Fi g. |q using (7.4) 
The adhesion energy is obtained from the definition (2.12]). Equation (4.5) with (5.3) implies 



a a d 



d/2 



dx 



[h'^ + eh"^ + ^-\h~h,,~2\x)^^ 



Using (7.4), the above integral can be calculated analytically yielding 



U ^Uq 4x2 fl, Iidv+,dr,^)+I{dr]^,dT]+) 



(7.7) 



(7.8) 



with 



I{u,v) 



2 ^(e/g.)M 4 ^[i + 2(e/e.) 

(^/d)2 |l+et (e/rf)2(w2_i;2) 

^^' ' ^ '^tanh(u/4) \ 



Fig. H (dashed lines) illustrates ( |7.8| ) and compares it with the exact numerical solution detailed below. 



(7.9) 
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B. The exact numerical solution 



In this subsection we pres ent a numerical solution which accounts for the full Van der Waals interaction (2.2) and 
the bending energy term in (2.11). One can functionally minimize (2.11) but the resulting Euler-Lagrange equation 
is non-linear and very complicated. Given that (7.2) implies the boundary conditions ( [7. 3D , the resulting problem is 
extremely awkw ard to tackle even numerically if one tries to solve the equation directly. Instead, we choose first to 
discretize ( 2.11 ) and then minimize the free energy functional with respect to all of the discrete variables. This has 
the advantage that the boundary conditions, (7.3), can be easily incorpor ated. 

For one dimensional structure, that is relevant to the trench geometry, (2.11) becomes 



T\h] = / dxl a\/lTh^ +- ■ — 
'1 2 (1 



//2 



- + V{h;z„0) 



(7.10) 



Due to the symmetry of (7.2), we only need consider a solution for < a; < d/2 and adopt a standard discretization 
process by dividing the interval, i.e. 



h'{xk) 



h"{xk) 



A2 



(7.11) 



where Xk = kA, k = 1, . . . , N and A = ■^. Here, N is the number of points makin g up t he o ne-dimensional lattice 
and typically was chosen somewhere between 100 and 400. For the simple surface (7.2), (2.3) can be broken into a 
(convergent) infinite sum of integrals each of which can then be evaluated analytically. 

Carrying out this procedure we find Figs, g and 0. From the adhesion energies plotted in Fig. ^, there is a region of 
good agreement between the Deryagin and exact solution. This agreement occurs for low d and A where (pj) holds. 
As the Deryagin approximation underestimates the attractive potential in this case, the adhesion energy is always 
less than the exact result and it is also reassuring that the exact U does not vanish for high d or A. This prediction 
of the Deryagin approach is clearly an artifact of going beyond the limits of perturbation theory. Generally, we see 
that the adhesion energy decreases with greater values of the roughness, i.e. large A or d (for small d any structure is 
effectively 'washed out' and not seen by the membrane). 

The degree of penetration of the membrane into the 'V'-shaped trenches is shown in Fig. u\ where the membrane 
height above the middle of the trench, Sh{Q) = h{0) — ho, is plotted. One can see that the membrane always lies 
further away from the substrate than if the latter were entirely planar and so there is no penetration into the surface 
indentations. However, this could be encouraged by having flat regions separating each trench (see Ref. M for a 
similar example). As the surface becomes rougher, it also becomes more repulsive {U decreases) and so the membrane 
moves outwards. It is perhaps surprising to see that the exact solution lies furthest away from the substrate d espite 
having a higher adhesion energy. This is likely to be a consequence of higher order bending energy terms in ( 7.10| ) 
reducing the amplitude of the membrane configuration and so increasing its height at the center of a trench. 



VIII. CONCLUSION 



A significant experimental question is whether or not substrate structure encourages membrane adhesion. An im- 
portant conclusio n of our study is that chemical structure always increases a substra te's attractiveness; the membrane 
potential energy (6.2) is clearly negative and also greater than the bendi ng e nergy, (6.3). Consequently, the adhesion 
energy increases. Rough surfaces are unfortunately more ambiguous as (3.8) is not of a definite sign. However, if the 
Deryagin approximation is invoked {K{q) +V ~ 1), then the potential energy contribution is always positive. There- 
fore, we would expect roughness to usually decrease a substrate's attractiveness and lead to a drop in the adhesion 
energy. We should emphasize that for surfaces for which non-local effects are important this may no longer be the 
case. 

Finally, some comments are pertinent on the validity of our approach. Both analytical methods breakdown when 
the amplitude of the structure, be it geometrical or chemical, becomes large. This is to be expected as our analysis 
is fundamentally a perturbation method and can only be confidently followed when (6.1) holds. The linear response 



technique is an improvement over Deryagin and is particularly appropriate for rough surfaces where the additional 
non-local effects lead to an increase in the amplitude of an adhering membrane. For smoothly varying surfaces 
these effects can even lead to the surface becoming attractive — a result that is not predicted by the local Deryagin 
approximation. Unfortunately, it is difficult to identify the particular geometries for which non-local Van der Waals 
contributions are important but for those surfaces of biotechnological interest, i.e. with trenches or indentations etched 
into them, they do not seem to lead to radically different behavior. To conclude, if one wishes an analytical guide to 
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how a certain substrate structure will influence membrane adhesion and if that structure can be conveniently described 
in Fourier space then the linear response description is the method of choice. Failing this the Deryagin approximation 
is quick and easy to apply if only normally adequate for small amplitude effects. 
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K = 35T CT = 1.7x 10"^ Jm"2 5 = 38 A T = 4.1 x 10" 

Ao = 2.6 X 10"^^ J 6 = 0.93Jm"2 a"^ = 2.2 A 

a ~ 49.3 A ho^ 0.61a ~ 30 A C/o ~ 0.298cr v 

^^ ~ 0.21a ^^ ~ 1.97a C - 18.62a V ^ 1.30 

TABLE I. The various parameters, chosen and calculated, for a supported membrane. For definitions see text. 
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FIGURE CAPTIONS 

Figure ffl Supported membranes on structured surfaces. In (a), the membrane is adhering to a rough but otherwise 
homogeneous surface. The reference p-plane is shown as a dashed hne. The height of the lower membrane Hpid 
leaflet and the surface, measured from this plane, are denoted by h{p) and Zs{p), respectively. Layered and 
columnar chemical structure are sketched in (b) and (c), respectively. 

Figure |^ A plot of the various interactions described in Sec. Ill, with parameter values given in Table g. Here all potentials 
are measured in units of the membrane tension, a = 1.7 x lO^^Jm^ , and lengths in terms of a ~ 49. 3A. The 
total potential has a minimum at ho ~ 0.6a. 

Figure M Membrane adhesion energy for a substrate which (a) is flat and has a sinusoidally varying Hamaker constant 
(Ac — 10.0) and (b) has a sinusoidal surface conflguration (Ag = 2.0a) plotted against wavenumber. The 
Deryagin solution is also shown with a dashed line. 

Figure Q A membrane adhering to (a) a substrate with a square wave surface conflguration and (b) to a flat substrate 
periodically striped with two different chemical compounds. 

Figure M A typical membrane configuration predicted by the Deryagin approximation for adsorption above a homogeneous 
substrate sculptured with 'V-shaped trenches. The parameters for a single trench are d = 10a and A = 0.05, 
which implies a width of around 500 A and depth of approximately 25 A. The adhesion energy is C/ ~ OASUq. 

Figure ^ A comparison of the Deryagin predictions (dashed lines) and an exact numerical solution (heavy lines) for the 
adhesion energy, U, above a substrate patterned with 'V shaped trenches. In (a) U/Uq is plotted as a function 
of d/a with A = 0.1, while for (b), it is shown as a function of A and d = 2a. For small d/a or A the Deryagin 
method provides a good approximation to the numerical result. 

Figure |^ The membrane height in the center of a 'V shaped trench measured with respect to its planar height, Sh{0) — 
h{0) — Hq. In (a) Sh{0) is plotted as a function of the trench width d/a while in (b) it is plotted as a function of 
the trench amplitude with d = 2a. The numerical solution (heavy line) soon departs from the Deryagin result 
(dashed line) with good agreement again only occurring for values of the adhesion energy close to Uq. 
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